Quantum radiation reaction force on a one-dimensional cavity with two relativistic 

moving mirrors 
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We consider a real massless scalar field inside a cavity with two moving mirrors in a two- 
dimensional spacetime, satisfying Dirichlet boundary condition at the instantaneous position of 
the boundaries, for arbitrary and relativistic laws of motion. Considering vacuum as the initial field 
state, we obtain formulas for the exact value of the energy density of the field and the quantum force 
acting on the boundaries, which extend results found in previous papers For the particular 

cases of a cavity with just one moving boundary, non-relativistic velocities, or in the limit of infinity 
length of the cavity (a single mirror), our results coincide with those found in the literature. 
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The Dynamical Casimir effect has been investigated 
since the 1970s 0-0], and has attracted growing atten- 
tion. It is related to problems like particle creation in 
cosmological models and radiation emitted by collaps- 
ing black holes d, M, decoherence |j[ , entanglement [l(| , 
the Unruh effect [ll|, among others. Models of a single 
mirror have been investigated and also cavities with one 
moving boundary have been studied in many papers (for 
a review see Ref. [HI). I n contrast, the problem of a cav- 
ity with two moving boundaries has been investigated re- 
cently and relatively few papers on this subject are found 
in the literature (for instance, Refs. [LML71 ]). A cavity 
with two oscillating mirrors can exhibit situations of con- 
structive and destructive interference in the number of 
created particles, depending on the relation among the 
phase difference of each boundary, the amplitudes and 
frequencies of oscillation [i~3l - [l5| . Ji, Jung and Soh [3] , 
considering the expansion of the quantizing field over a 
instantaneous basis and a perturbative approach, investi- 
gated the problem of interference in the particle creation 
for a one-dimensional cavity with two boundaries moving 
according to prescribed, non-relativistic and oscillatory 
(small amplitudes) laws of motion. Dalvit and Mazzitelli 
[i~5| extended the field solution obtained by Moore [|[ 
for the case of a one-dimensional cavity with two moving 
boundaries, deriving a set of generalized Moore's equa- 
tions, also obtaining the expected energy-momentum ten- 
sor for this model, generalizing the corresponding for- 
mula obtained by Fulling and Davies In Ref. 15 1 
the set of generalized Moore's equations was solved for 
the case of a resonant oscillatory movement with small 
amplitude, via renormalization-group procedure. Li and 
Li [2| applied the geometrical method, proposed by Cole 
and Schieve [Hj], to solve exactly the generalized Moore 
equations obtained by Dalvit and Mazzitelli [l5[ , and also 
used numerical methods to obtain the behavior of the en- 
ergy density in a cavity for particular sinusoidal laws of 
motion, with small amplitude [3] . On the other hand, as 
far as we know, there is no paper in literature devoted 
to obtain formulas which enable us to get directly ex- 



act values for the quantum force and energy density in 
a nonstatic cavity for arbitrary laws of motion for the 
moving boundaries, including non-oscillating movements 
with large amplitudes, which are out of reach of the per- 
turbative approaches found in the literature. 

In the present paper we consider a real massless scalar 
field satisfying the Klein-Gordon equation (we assume 
throughout this paper h = c = 1): (<9f — <9 2 ) <p {t, x) = 0, 
and obeying Dirichlet conditions imposed at the left 
boundary located at x = L(t), and also at the right 
boundary located at x = R(t), where L(t) and R(t) are 
arbitrary prescribed laws of motion, with R(t < 0) = Lq 
and Lit < 0) = 0, where Lq is the length of the cavity 
in the static situation. Considering vacuum as the initial 
field state, we obtain formulas for the exact value of the 
energy density of the field and the quantum force acting 
on the boundaries, showing that the energy density in 
a given point of the spacetime can be obtained by trac- 
ing back a sequence of null lines, connecting the value of 
the energy density at the given spacetime point to a cer- 
tain known value of the energy at a point in the "static 
zone", where the initial field modes are not affected by 
the disturbance caused by the movement of the bound- 
aries. Our formulas generalize those found in literature 
[![, where this problem is approached for a cavity with 
only one moving mirror. For the particular cases of a 
cavity with just one moving boundary, non-relativistic 
velocities, or in the limit of large length of the cavity (a 
single mirror), our results coincide with those found in 
the literature [H 11 ll~9l| , 

Let us start considering the field operator, solution of 
the wave equation, given by [1 51 ] 
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where the field modes are 



(t, x) 



(i) 
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with v = t + x, u = t — x, and 

G[t + L (<)] — F[t — L (*)] = (2a) 
G [t + J? (t)] - F [t - R (£)] = 2. (2b) 

The set of equations fl2|) , obtained by Dalvit and Mazz- 
itelli exploiting the conformal invariance of the model 
15], is a generalization of the Moore equation Q, which 
can be recovered doing L(t) = in these equations. The 
renormalized energy density in the cavity is given by [l5j ] 



where 
fc(z) 
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In the present paper we use Eqs. (|2a|) . (|2bj) . (|4a[) and 
(I4b[) to obtain the following set of equations for the func- 
tions fa and /i?: 



f G [t+R(t)} = f F [t~R(t)]A R (t)+B R (t), 
f G [t + L(t)] = f F [t-L(t)]A L (t)+B L (t), 



with 



A q (t) 



1 + 9' (*) 



(5a) 
(5b) 
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B g (t) = - 
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q'» (t) 
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?" 2 (*)</' (t) 



47T [1 + g/ (t)f[l- 9 ' (t)]< 



(7) 



where, hereafter, q can represent R or L. Eqs. (|5a[) and 
(|5bj) are an extension of the corresponding equation for 
/, valid for a cavity with just one moving mirror, found 
in Ref. 0]. If we consider the particular case of L(t) = 
in Eq. ([5]), we recover the corresponding result found in 
Ref. 0| . For (t < 0) we have f G (v) = f F (u) = f (s) = 
7r/(48Lg), and (Too) = - 7r /(24Lo): which is the Casimir 
energy density for this model. 

Now, our aim is to solve the Eqs. (f5a| and (|5b|) recur- 
sively, using a geometrical point of view. Let us exam- 
ine the cavity in the nonstatic situation (t > 0). The 
field modes in Eq. ([T]) are formed by left and right- 
propagating parts. As causality requires, the field in 
region I (v < L a and u < 0) (see Fig. [1]) is not af- 
fected by the boundaries motion, so that, in this sense, 




FIG. 1: Boundaries trajectories (solid lines) . The dashed lines 
are null lines separating region I from II and III, and these 
ones from region IV. 



this region is considered as a "static zone" . In region II 
(v > Lq and u < 0), the right-propagating parts of the 
field modes remain unaffected by the boundaries motion, 
so that region II is also a static zone for these modes. 
On the other hand, the left-propagating parts in region 
II are, in general, affected by the boundary movement. 
Similarly, in region III (it > and v < Lq), the left- 
propagating parts of the field modes are not affected by 
the boundaries motion, but the right-propagating parts 
are. In region IV (v > Lq and u > 0), both the left and 
right-propagating parts are affected. In summary, the 
functions corresponding to the left and right-propagating 
parts of the field modes are considered in the static zone 
if their arguments are, respectively v < Lq and u < 0. 
Then, we have fc(v < L ) = f (s 1 and f G (u < 0) = f^ s \ 

For a certain spacetime point ft, x) , the energy tensor 
( Too is known if its left and right-propagating 

parts, taken over, respectively, the null lines v = Z\ and 
u = zi (where Z\ = t + x and Z2 = t — x), are known; or, 
in other words, (Too (^^)) is known if f G (v)\ v _ Zl and 
f F (u)\„ are known. Li and Li Q used a recursive 
method [l8| to obtain the functions G and F for 
general laws of motion of the boundaries, tracing back a 
sequence of null lines until a null line gets into the static 
zone where the G or F functions are known. Here, we 
adopt this method to obtain f G and f F , extending the 
work done in Ref. . Let us assume that (t, x) belongs 
to region IV, and that the null line v = z\ intersects 
the moving mirror trajectory at the point [ti,R(ti)] 
(see Fig. |2(a)[ ), so that t + x = t\ + R(ti). We have 
R(h)]. Using the Eq. ([Sa|). we 



fc(v)\ v=zi = f G [h 
get f G [ii + R (ti)] = f F [h - R (t x )] A R (t x ) + B R (t x ). 
If ti - R(ti) < 0, then the null line u_=J i - R(h) 
is already in the static zone (Fig 
we can write f F [ti — R (ti)] = 

f G [ti + R(h)} 



./ 



, so 
and 



that 
also 



f^A R (h) + Buih), and we 
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FIG. 2: Sequence of null lines (dotted lines) connecting a 
point (t, x) to a static zone. The dashed lines are null lines 
separating region I from II and III, and these ones from region 
IV, as presented in Fig. [U In Fig. |2(a)l we see the case of 
one reflection (jig = 1), whereas in Fig. 2(b) we see the case 
n G = 2. 



can say that the number of reflections uq to get into 
the static zone is, in this case, rig = 1. On the 
other hand, if t\ — R(t-y) > (case shown in Fig. 
|2(b)[ ) wc can draw another null line v = t 2 + L (t 2 ) 
intersecting the world line of the left boundary at 
the point [t 2 ,L(t 2 )], with t x - R(t x ) = t 2 - L(t 2 ). 
In this case we have, using (|5b"f . fa [tx + R (ti)] = 
{fa [*2 + L {t 2 )\ - B L (t 2 )} A R (t x ) /A L (t 2 ) _^Br (f x ). 
If 



*2 

/G ]*2 - 
{/ (S) - 

if t 2 ^ 

v = t 2 - 



L{t 2 ) < L (see Fig. 



2(b) I, then 



i(«a)] = f {s \ 

B L (t 2 )} A R (ti) /A L (i a ) + Br (tx) and n G = 2. 
■ L (fa) > ^Oj we assume that the null line 
(-£ (fa) intersects the right boundary at the point 
[i 3 ,i?(t 3 )], then t 2 + L(t 2 ) = t 3 + R(t 3 ) and we get 
u = £3 — R{t 3 ). We repeat this procedure up to a null 
line gets into a static zone, where the function fp or /g 
is known. In summary, we obtain for Jq: 



f G (z) = /«I G (z) + B G (z) , 



where, for tiq (z) even, we have 



A, 



■w= n 



k=0 



(1 - 4,0) 



Ar [*2fc-l(-z)] 



3fc,0 
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(9a) 
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^R [*2fc-l(z)] [*2fc(z)] 

AR[t 2fc _i(z)] 



-B L [t 2k (z)}j 



n 



Ar[%_i(^)] 



(9b) 



with S symbolizing Kronecker's delta function. For 
tiq (z) odd we have 



A G (z) 



Bg (z) 



*g(*)-i 

2 

n 



Ar [*ak+i(*)] 



(l-4,o)^L [<2fe(^)]+4,0 



, (10a) 



2 

E [fek+i(*)] - (!- 4o) £?l M*)]] 



fc=0 

fc 

><n 

3=0 



, v ^rN-i(^)] x 



(10b) 



Note that the number hq of reflections and the sequence 
of instants t±, t nG depend on the argument z. The set 
of instants mentioned in Eqs. © and (|10p are calculated 
via [2j: 



*2i+l — R(t2l+l) 
t 2 l+2 + L{t 2 i +2 ) 



h + R(h), 

t 2 l+2 — L(t 2 i +2 ), 
t 2 l+3 + R{hl+3), 

0,1,2... 



(11) 



To solve recursively the set of equations §5§ for fp, 
we start assuming that the null line u = t — x in- 
tersects the worldline of the left mirror at the point 
[t~i, £ (f~l)], so that t — x = i\ — £(t~i). Thus we have 
fp ( u )L=z 2 _= If [tx - L (ti)]. Using the Eq. (50, we get 
/f [fx ~ £ (tx)] = {/g [tx + £ (tx)] - B L (ti)} /Ai (ti). If 
<i + L(t\) < Lq, then the null line v = t\ + L(ti) 
is already in the static zone, so that we can write 
fclh+Lih)} = and also fp [i\ — L (ti)] = 

|j( s ) _ b l (t\)} /Al (tx), and we can say that the num- 
ber of reflections rip to get into the static zone is, in this 
case, np — 1 (see Fig. 3(a)). On the other hand, if 
ti + L (ti) > Lo (as shown in Fig. |3(b)[ ) we need to find 
fa [tx + L (ti)] recursively via Eq. In general, we get 



where 



f F (z) = f^A F (z) + B F (z), 
A G {t 1 {z) + L[t 1 {z)}} 



A F (z) 



AlUz)] 



(12) 



(13a) 



^^Mw+pi^i, ( 13b ) 



AlUz)] 

with the function t\{z) calculated via 
z = ti-L(ti). 



(14) 



The formulas (|9]), (| 10[) and ([T3|) generalize those for A 
and £? found in Ref. [1], which are valid for a cavity with 
just the right boundary in movement. 
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FIG. 3: Sequence of null lines (dotted lines) connecting a 
point (t, x) to a static zone. The dashed lines are null lines 
separating region I from II and III, and these ones from region 
IV, as presented in Fig. [TJ In Fig. [3(a)] we see the case of 
one reflection (uf — 1), whereas in Fig. 3(b) we see the case 
n F = 2. 



From Eqs. ©, ^ and (fT2"1) . we get the exact formula 
for the renormalized energy density as 



<T 00 (t.a;)) --/ (s) A G (v)+A F (u) 



-B G (v)-B F (u). 

(15) 

Eq. (fT5")) gives directly the exact values for the energy 
density in a nonstatic cavity for arbitrary laws of motion 
R(t) and L(t). Since Too = T\\ in this model, we have the 
following exact formulas for the renormalized quantum 
forces J-g = (Too [*, #(*)]) and T L = - (T 00 [t, L(t))) (see 
Refs. [13, acting, respectively, on the right and left 
boundaries: 

F R (t) = -f {s) [A G [t + R(t)]+A F [t-R(t)]} + 

-B G [t + R(t)] — Bp [t — R(t)} , (16) 

•Fi(t) = / (s) {i G [t + T(t)]+i F [i-T(t)]} + 

+£ G [t + L(t)] + Bf. [f - L(t)] . (17) 

Next we examine the behavior of these forces in each 
region pointed in Fig. [TJ 

In region I (Fig. [TJ), we have n G = np = 0. Then, 
Eqs. © and flT3]) give: A G (z) = Ap (z) = 1 and 
B G (z) = B F (z) = 0. This results, as expected, in the 
static Casimir force 



[Cas] 



= = -*/(24L8), 



acting on the boundaries. 

In region II, we have ng = 1 and n F = 0. For this 
case, Eq. ([T3")) gives (w) = 1 and Bp (u) — 0, whereas 
from Eq. (fTTjj) we have A G (v) = and B G (v) = 

Br [ti(v)]. To calculate the force in Eq. (|16l) we do 



v — s- t + -R(i), and obtain ti(i>) as already discussed: t + 
R(t) = ti+R(h) =>h=t. Then we get A G [t + R (t)] = 
A fl (i) and B G [t + R (i)] = 5;?. (i). The force J" fl (t) on 
the right boundary in region II, now relabeled as J~ R ^ (t) 
is 



^i II) W--/ ls, [1 + Ar (t)]-B R (t). 



(18) 



From this formula, we can obtain an analytical result for 
an arbitrary law of motion R(t). Note that in Eq. (|T8j) 
the subscript T is not found, since the quantum force for 
the worldline in region II has no influence of the move- 
ment of the left boundary. Considering the limit To ~~ ^ o° 
we recover the quantum radiation force !Fq correspond- 
ing to the unbounded field, acting on the left side of a 



single mirror: limi, _>. oc , F R l \t) — J 7 ^' (t), where 



7W1 



A~ u Ht) = -B q (t). 



(19) 



we get ^(t) 



In the non-relativistic limit, from (fP? 

Fr 0,8 ^ + i?/(127r), and adding tne limit l,q — > 00 
we recover the approximate quantum radiation force 
T R 11 ^ (t) ss R/ (127t), which acts on the left side of a single 
mirror [2lj . 

In region III, we have n G — and np = 1. For 
this case, Eqs. © and (US]) give Ap (u) = 1/A L [ii (u)]; 
Bp (u) = -B L [h (u)J /A L [tj (u)]; i G (w) =1; £ G (u) = 
0. Considering u — > t — L(t) and t — L(t) =t\ — L(t\) => 
<i = t, the force Jz,(t) on the left boundary in this region, 
now relabeled as J-^ 11 ^ (t) is 



Considering the limit T 



1 



A L (t) J A L (t) ■ 



(20) 



00 we recover the quan- 
tum radiation force Tq corresponding to the un- 
bounded field, acting on the right side of a single mirror: 

limzo-voo ^(t) = F { L +U) (t), where 



gg (*) 

A q (*)' 



(21) 



From Eqs. (|T9j) and (|21j) we recover the total quantum 
force Fq (t) acting on a single mirror at vacuum, with a 
prescribed trajectory x = q(t): 



= (l + f){[?q/(2ir)]/(l-q>y 
+ [9/(6^)]/(l-9 2 ) 3 }, 



which is in agreement with that found in literature (see 
Ref. [Hj]). In the non-relativistic limit, we reobtain the 
approximate quantum radiation force Jq\t) ~ <7 /(67r) 



To compute the total forces !F R tot ^ and j:^ ULI acting 
on, respectively, the right and left boundaries, for any of 



-(tot) 
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the regions II, III or IV showed in Fig. [TJ we need, in 
addition to Eqs. (TH)1) and (JTTJ), to take into account the 
remaining dynamical Casimir forces corresponding to the 
vacuum field outside the cavity, which are given by Eqs. 
(fHI and (Ull)- We write: 



£*> = ^(t)+^(t), 



(22) 
(23) 



Eqs. (j2"2")l and (f2"3")l enable us to calculate directly and an- 
alytically the total quantum forces acting on both mirrors 
for arbitrary laws of motion R{t) and L{t), in regions II 
or III, because for these regions Eqs. ([To]) and (IT71) are 
replaced by their particular cases given by Eqs. (I18|) and 

In region IV (see Fig. [T]), in general it is difficult to 
obtain exact analytical results for the quantum forces 
(HHJ) and (|17[) . for arbitrary trajectories R(t) and Lit). 
The difficulty is in solving equations like t\ — R(t\) — 
t,2 — L (£2) (see Eq. (jlip h which arise after a second re- 
flection (no > 2 or/and > 2). Trajectories can be 
constructed to give analytical solutions to these equa- 
tions, but a large class of relevant laws of motion do not 
result in exact analytical solutions. However, our results 
enable us to obtain exact numerical results for the quan- 
tum force acting on the moving boundaries of a cavity for 
an arbitrary law of movement, including non-oscillating 
movements with large amplitudes, which are out of reach 
of the perturbative approaches found in the literature, as 
we will examine next. In this context, let us apply our 
formulas to the following particular non-trivial trajec- 
tory, which is based on the one proposed by Haro in Ref. 
221: 



L{t) = kl ln[cosh(t)], 
R(t) = Lq + Kflln[cosh(t)]. 



(24a) 
(24b) 



Considering, for instance, kr = —kl = 0.1 (Fig. 4(a) I, 



FIG. 5: The solid line shows the total force J-r (t), whereas 
the dashed line shows the force J-n{t), both for the law of 
movement (|24[) . with kr — — kl = 0.1 and Lq = 1. 



|4(b)[ ), we have the mirrors in movement with relativistic 
velocities, but keeping constant the cavity length. 

In Figs. [5] and [I using our formulas (|5])-([T3"]l and ([2"2]l. 
we plot the time evolution of the quantum force (t) 
and Tr (£ ), for, respectively, the cases kr — —Kl = 0. 1 
(see Fig 



4(a) I, and kr = kl = 0.1 (see Fig. 



we have an expanding cavity with large amplitude and 
relativistic velocities. If we consider kr = kl = 0.1 (Fig. 



We can see discontinuities of the derivatives for J- R tot ^ 
and Fr (t) . These discontinuities always occur when 
the front of the wave in the energy density meets the 
right boundary. In the case, for instance, showed in 
Fig. [5j when t — the left boundary starts to move 
and generate a wave in the energy density, propagating 
rightward and meeting the right boundary at the instant 
t = Ti PS 1.05, calculated via equation t\—R (ri) = 0, and 
which corresponds to the first discontinuity of the deriva- 
tive showed in Fig. [SJ At t = 0, another front of wave 
is generated by the right boundary, propagating leftward 
and meeting the left boundary at the instant t\ ps 1.05, 
and then reflected back and meeting the right boundary 
at the instant T2 as 2.25, calculated from the equation 
T2 — R (T2) = Ti — L (ri). This instant corresponds to the 
second discontinuity of the derivative showed in Fig. [S] 
Since the length of the cavity remains the same in the 
case showed in Fig. |4(b)| the quantum force J 7 ^* *^) os- 
cillates around the static Casimir force (Fig. [S]) , whereas 
it goes to zero in the case showed in Fig. (Fig. [5]) where 
the boundaries go to an asymptotic behavior of infinity 
lenght and constant velocity. 

Summarizing our results, the formulas obtained in the 
present paper enable us to get directly exact values of the 
energy density of the field and the quantum force acting 
on the boundaries of a nonstatic cavity for arbitrary laws 
of motion for the moving boundaries, for vacuum as the 
initial state of the field. Eqs. (|5al) and (I5b|) are an exten- 
sion of the corresponding equation for a cavity with just 
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-0.115- /\ 




FIG. 6: The solid line shows the total force T R tot \t), the 
dashed line shows the force Tr^), both for the law of move- 
ment (|24[) . with he = hl =0.1 and Lq = 1. The dotted line 
shows the static Casimir force . 



one moving boundary found in Ref . , and the achieve- 
ment of fc and fp recursively, tracing back a sequence 
of null lines, can be viewed as an extension of the work 
done in Ref. @. Formulas ©, (flU]) and (P|) general- 
ize those found in Ref. For the particular cases of 
a cavity with just one moving boundary, non-relativistic 
velocities, or in the limit of infinity length of the cavity 
(a single mirror), our results are in agreement with those 
found in the literature [J II H, d ■ 

The present results enable investigation of several 
problems (usually treated by perturbative approaches in 
the literature) with an exact approach and also out of the 
regime of small amplitudes. For instance, those related to 
the inertial forces in the Casimir effect with two moving 
mirrors [l7| . or the interference phenomena in the pho- 
ton production [l4j • These issues are under investigation 
and will be discussed in future papers. 
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